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Samenvatting 

Een  veel  gebruikt  veiligheidscriterium  voor  de  opslag  van  kruit  is  de  veilige  diameter.  Een  opgesla- 
gen  pakket  kruit  mag  deze  diameter  niet  overs chrijden,  opdat  de  warmte  die  door  langzame  ont- 
ledingsreacties  wordt  geproduceerd,  snel  genoeg  naar  de  omgeving  afgevoerd  kan  worden.  Dit 
criterium  wordt  gewoonlijk  toegepast  op  de  afineting  van  een  enkele  verpakkingseenheid.  De  af- 
koeling  kan  echter  belemmerd  worden  door  een  dichte  stapeling  van  verpakte  eenheden  kruit  in 
een  opslaggebied,  waardoor  het  systeem  toch  instabiel  kan  worden. 

Om  dit  te  onderzoeken  is  een  stapel  van  twee  lagen  van  vijfentwintig  dicht  op  elkaar  gestapelde 
vaten  kruit  gemodelleerd  met  een  eindige  elementen  methode  en  de  resultaten  zijn  vergeleken  met 
het  gedrag  van  een  enkel  vat.  De  warmtestroming  in  de  ruimte  tussen  de  vaten  is  gemodelleerd 
met  behulp  van  een  eftectieve  warmtegeleidingsco££Sd£nt  van  1  W/mK. 

De  stapel  had  een  10  tot  15*C  lagere  kritische  omgevingstemperatuur  dan  een  enkel  vat.  De  laag- 
ste  kritische  omgevingstemperatuur  die  is  gevonden  was  circa  50’C.  De  (effectieve)  warmtegelei- 
dingscoSfiSdfinten  van  het  kruit  en  de  lucht  tussen  de  vaten  zijn  de  belangrijkste  parameters  die  de 
thermische  stabiliteit  van  de  stapel  bepalen. 

Summary 

A  commonly  used  safety  criterion  for  the  storage  of  gunpowder  is  the  critical  diameter.  This  dia¬ 
meter  must  not  be  exceeded  by  a  stored  package  of  gunpowder  to  ensure  that  the  heat  production 
caused  by  slow  decomposition  reactions  can  be  compensated  by  loss  of  heat  to  the  surroundings. 
This  safety  criterion  is  normally  applied  to  the  size  of  a  single  package.  However,  the  cooling  effect 
can  be  decreased  by  the  close  stacking  of  packages  in  a  storage  area,  which  could  render  die 
system  unstable. 

A  stack  of  two  layers  of  twenty-live  closely  packed  drums  of  powder  was  modelled  using  a  finite 
element  method  and  the  results  were  compared  with  the  behaviour  of  a  single  drum.  The  heat  flow 
in  the  space  between  the  drums  was  modelled  using  an  effective  conductivity  of  1  W/mK. 

The  stack  had  a  10  to  15*C  lower  critical  ambient  temperature  as  compared  to  the  single  drum. 
Critical  ambient  temperatures  as  low  as  approximately  50*C  were  encountered.  The  (effective) 
conductivity  of  the  powder  and  the  air  between  the  drums  are  the  most  important  factors  that 
determine  the  thermal  stability  of  the  stack. 
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1  INLEIDING 

Bij  de  normale  opslagtemperatuur  venoont  kruit  een  zeer  langzame,  maar  met  verwaarloosbare, 
ontiedingsreaccie.  Hie  rbij  komt  warm  re  vrij,  waardoor  het  kruir  in  temperatuur  stijgt.  Als  deze 
warmte  niet  snel  genoeg  wordt  afgevoerd,  is  het  gevaar  van  spontane  ontbranding  aanwezig.  De 
omvang  van  een  pakket  kruit  in  de  vorm  waarin  het  wordt  opgeslagen,  bepaalt  mede  of  de  gepro- 
duceerde  warmte  goed  kan  worden  afgevoerd. 

De  stabiliteit  van  kruit  kan  men  bepalen  door  de  warmteontwikkeling  van  het  kruit  te  meten  [1]. 
Daarmee  wordt  de  veilige  diameter  bepaald,  die  niet  overschreden  mag  worden  voor  een  opgesla¬ 
gen  eenheid  kruit.  In  de  prakdjk  liggen  deze  eenheden  gestapeld  in  een  bepaalde  verpakking. 
Hoewel  de  eenheden  zelf  een  veilige  diameter  hebben  is  het  de  vraag  in  hoeverre  dat  criterium  nog 
van  waarde  is  bij  stapeling.  De  buitenste  eenheden  vormen  namelijk  een  soort  "warmte-isolatie" 
voor  de  binnenste  eenheden. 

In  opdracht  van  Bureau  Hinderwetzaken  van  de  Direcde  Materieel  van  de  Koninklijke  Landmacht 
is  de  invloed  van  de  stapeling  van  kruit  theoredsch  onderzocht.  Met  behulp  van  het  eindige  ele- 
menten  programma  "Abaqus"  is  een  computermodel  opgesteld  van  cilindervormige  vaten 
(fiberdrums)  die  in  twee  lagen  van  vijf  bij  vijf  vaten  zijn  opgestapeld.  Een  be  perking  van  het  model 
is,  dat  de  stroming  van  lucht  tussen  de  vaten  niet  uitgerekend  kan  worden.  Daarom  is  het  warmte- 
transport  dat  door  luchtstromingen  veroorzaakt  wordt  benaderd  door  het  te  beschrijven  als  een 
warmtegeleidingsmechanisme . 

Door  middel  van  simulades  is  nagegaan  hoe  hoog  de  omgevingstemperatuur  kan  zijn,  zonder  dat 
er  een  wegloopreactie  optreedt.  Dit  is  vergeleken  met  de  situaue  van  een  los,  niet-gestapeld  vat 
kruit.  Daarbij  is  gekeken  naar  het  temperatuurverloop  in  de  tijd,  en  naar  het  temperatuurprofiel 
dat  uiteindelijk  in  de  stapel  omstaat,  als  de  sitvatie  thermisch  stabicl  is. 


In  de  eerste  stap  is  de  stapel  met  een  eenvoudig,  tweedimensionaal  model  beschreven.  In  de 
tweede  stap  is  het  model  uitgebreid  tot  een  volledig  driedimensionaal  model.  Tevens  is  de  gevoe- 
ligheid  van  het  model  voor  de  waarde  van  de  gebruikte  fysische  parameters  onderzocht. 
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2  THEORIE 

2.1  Thermlsche-explosietheorie 

In  het  hier  gebruikte  model  wordt  uitgegaan  van  een  reactiekinetiek  volgens  vergelijking  (1).  Dit  is 
een  Arriienius  vergelijking  voor  een  nulde-orde  reactie. 


q  -  Qo  exp(-E/RT)  [W/kg]  (1) 

Hierin  is  q  de  warmteproduktie  [W/kg],  Qo  de  warmteproduktiefactor  [W/kg],  E  de  active- 
ringsenergie  [J/mol],  R  de  gasconstante  Q/mol  K]  en  T  de  absolute  tempera tuur  [K]. 

Het  warmtetranspon  door  een  materiaal  wordt  bepaald  door  de  wanntegeleidingsco€fficiCnt  en  de 
tempera  turn-gradient  in  het  materiaal. 


j=_lSw/m2]  (2) 

J  is  de  warmteflux  [W/m2],  1  is  de  warmtegeleidingscoCffidCnt  [W/mK]  en  dT/dx  de  gradient  van 
de  temperaruur  [K/m]. 

De  warmteoverdracht  van  het  kruitvat  naar  de  lucht  er  omheen  wordt  beheerst  door  vergelijking 
(3). 


J  =  h(Tk-TA)  [W/m2]  (3) 

Hierin  is  h  de  warmteoverdrachtscoefficient  [W/m2K]  en  Tk  en  TA  zijn  respectievelijk  de  tem- 
peratuur  van  de  buitenkant  van  het  kruitvat  en  de  tempera  tuur  in  de  directe  omgeving  van  het  vat. 

Het  verloop  van  de  temperatuur  in  de  tijd  volgt  uit  de  energiebalans: 


4r  =  a  £-y  +  Cp-q[K/s] 


a 


(4) 
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Voor  het  een-dimensionale  gevai  leidt  deze  tot  de  volgende  particle  differentiaalvergelijking.  t  is  de 
tijd  [sj  en  x  is  de  ruimtelijke  coOrdinaat  [m],  De  thermische  diffusiviteit  a  is  gedefinieerd  als 
a  =  l/(p.Cp);  p  is  de  dichtheid  [kg/m3],  Cp  is  de  soortelijke  warmte  [J/kgK]  en  q  is  de 
warmteproduktie  [W/kg], 

Met  de  Thermische  Explosie  Theorie,  zoals  beschreven  in  [2],  kan  aangetoond  worden  dat  een 
dergelijk  systeem  instabiel  is  als  de  omgevingstemperatuur  boven  een  bepaalde  kridsche  waarde 
ligt.  Voor  een  aantai  eenvoudige  geometrieCn  kan  de  kridsche  omgevingstemperatuur  uitgerekend 
worden.  Voor  een  oneindig  lange  cilinder  gaat  dat  als  volgt.  De  dimensieloze  grootheid  d  wordt 
gedefinieerd  met  veigelijking  (5)  (r  is  de  straal  [m]). 


d  =  (E/R  Ta2)  (r2  p/1)  exp(-E/R  TA)  (5) 

Volgens  de  Thermische  Explosie  Theorie  is  een  systeem  stabiel,  indien  d  beneden  de  kridsche 
waarde  (do)  ligt.  In  figuur  2.5  van  Bowes  [2]  wordt  dc  gegeven  als  functie  van  het  Biot-getal 
(hr/1).  Hiermee  kan  de  kridsche  omgevingstemperatuur  berekend  worden. 

2.2  Warmtetransport  in  een  vat  met  kruit 

Voor  een  ingewikkelde  geometrie,  zoals  een  stapel  vaten,  bestaat  er  geen  formule  om  de  kridsche 
omgevingstemperatuur  uit  te  rekenen.  Bovendien  is  het  warmtetransport  mssen  de  vaten  kruit  veel 
complexer  dan  het  warmtetransport  door  het  kruit  zelf.  In  de  tussenruimtes  tussen  de  vaten  treedt 
niet  alleen  warmtegeleiding  op,  maar  ook  convecue.  De  lucht  vertoont  een  complexe  stroming  en 
zorgt  daarmee  voor  warmte-uitwisseling  tussen  hete  en  koelere  plaatsen.  Deze  luchtstromingen 
worden  mede  veroorzaakt  door  de  dichtheidsverschillen  die  het  gevolg  zijn  van  temperatuur- 
verschillen.  De  mate  van  warmte-uitwisseling  tussen  de  vaten  is  daarom  niet  op  voorhand  bekend, 
maar  zal  liggen  tussen  twee  uitersten: 

-  de  verpakkingen  (en  tussenruimtes)  zijn  zeer  sterk  warmtegeleidend.  In  dat  gevai  zal  de  stapel 
geen  noemenswaardige  warmteweerstand  vertonen,  zodat  de  verpakkingen  met  kruit  als  losse 
eenheden  beschouwd  kunnen  worden,  zonder  invloed  van  de  stapeling; 

-  de  verpakkingen  (en  tussenruimtes)  zijn  sterk  isolerend.  In  dat  gevai  kan  het  kruit  zijn  warmte 
niet  kwijt,  zodat  de  temperatuur  tot  een  gevaarlijk  niveau  zal  oplopen,  ongeacht  de  grootte  van 
de  stapeling. 
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2.3  Eindige  elementen  methode 

De  stapel  wordt  gesimuleerd  met  behulp  van  het  softwane-pakket  Abaqus,  dat  gebaseerd  is  op  de 
eindige  elementen  methode.  Varianten  van  deze  methode  zijn  toepasbaar  op  verschillende  soorten 
problemen,  bijvoorbeeld  sterkteberekeningen,  stromingsproblemen  en  warmtetransportproblemen. 
Het  principe  is,  dat  het  te  onderzoeken  systeem  wordt  opgedeeld  in  een  groot  aantal  eenvoudige 
elementen,  zoals  kubussen,  tetraSders  en  dergelijke.  Deze  elementen  belnvloeden  elkaar  via  hun 
gemeenschappelijke  grensvlakken  en  hoekpunten.  Samen  vormen  al  deze  elementen  een  reken- 
rooster.  Als  de  materiaaleigenschappen  en  de  begin-  en  randvoorwaarden  van  alle  elementen 
gedefinieerd  zijn,  kan  het  gedrag  van  het  systeem  gesimuleerd  worden.  In  ons  geval  wordt  uit- 
gerekend  hoe  op  elk  roosterpunt  de  tempera tuur  in  de  tijd  verloopt. 

3  BESCHRIJVING  VAN  HET  MODEL 


3. 1  Beschrijving  van  de  opslagsituatie 

In  een  opslagruimte  is  kruit  opgeslagen  in  fiberdrums.  Dit  zijn  cilindervormige  vaten  van  een 
kanonachtig  materiaal  met  een  houten  deksel  (6  mm  dik  "plywood"  =  multiplex)  dat  met  een 
klemband  bevestigd  is.  Het  kruit  heeft  een  binnenverpakking  van  textiel  met  een  plastic  zak.  Er 
bestaan  verschillende  formaten.  In  het  model  is  uitgegaan  van  de  maten  van  een  veel  gebruikt 
type: 

diameter  355  mm 

hoogte  630  mm 

wanddikte  2,3  mm 

Er  wordt  aangenomen  dat  het  vat  100%  gevuld  is.  Er  is  geen  aparte  binhen-  en  buitenverpakking 
gemodelleerd,  maar  de  verpakking  wordt  in  het  model  door  tin  warmteoverdrachtscoCfficiint 
gekenmerkt.  De  waarde  van  deze  warmteoverdrachtsco£ffici£nt  is  6.0  W/m^K. 

De  vaten  staan  gestapeld  op  pallets  die  echter  niet  in  het  model  meegenomen  zijn.  Het  model  gaat 
uit  van  een  constante  temperatuur  in  de  opslagruimte  (buiten  de  stapel).  Stralingseffecten  zijn  niet 
meegenomen,  omdat  de  ervaring  leert,  dat  die  bij  lage  tempera turen  een  uiterst  geringe  bijdrage 
aan  het  watmtetransport  leveren.  Het  model  gaat  uit  van  een  stapel  van  vijf  rijen  van  vijf  vaten,  die 
dicht  tegen  elkaar  aan  gezet  zijn.  Er  zijn  twee  lagen  boven  elkaar  gestapeld,  zonder  tussenruimte 
(zie  figuur  3). 
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3.2  Beschrijving  van  het  rekenrooster 

In  de  eindige  elementen  methode  wordt  het  te  onderzoeken  systeem  opgedeeld  in  een  aantal 
elementen  van  eindige  grootte.  Hierdoor  wordt  een  rooster  (of  "grid")  gevormd,  aan  de  hand 
waarvan  het  systeem  wordt  doorgerekend.  Daarbij  wordt  gebruik  gemaakt  van  de  symmetrie  in  het 
systeem,  zodat  slechts  een  representatief  deel  van  het  systeem  hoeft  te  worden  doorgerekend. 

Voor  dit  onderzoek  zijn  enkele  verschillende  grids  gebruikt. 

In  de  eerste  plaats  is  een  enkel  vat  gemodeUeerd  in  een  twee-dimensionaal  model,  waarbij  wordt 
aangenomen  dat  alle  warmteafvoer  in  horizontale  richdng  plaats  vindt.  Dit  komt  in  rekenkundige 
zin  overeen  met  de  warmtestroming  in  een  oneindig  lange  dlinder.  Er  zijn  hiervoor  drie  grids  ge¬ 
bruikt  met  een  verschillende  mate  van  verfijning  (zie  figuur  1): 

-  grove  structuur,  acht  elementen; 

-  verfijnde  structuur,  twaalf  elementen; 

-  extra  fijne  structuur,  twee£ndertig  elementen. 


iron 
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Figuur  2 


2-D  model  voor  vijf  bij  vijf  vaten 
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Figuur  3  3-D  model  voor  twee  bij  vijf  bij  vijf  vaten 

Aangenomen  dat  alle  warmteafvoer  langs  de  zijwanden  plaatsvindt.  Dit  komt  overeen  met  een 
oneindig  hoge  stapel. 

Vervolgens  is  van  dit  model  een  vereenvoudigde  versie  gemaakt,  door  aan  te  nemen  dat  de 
warmteoverdrachtsco€fifid£nt  tussen  vaten,  en  de  luchtlagen  tussen  de  vaten  oneindig  hoog  is. 
Daardoor  konden  de  elementen  voor  warmteoverdracht  geSlimineerd  worden.  Tenslotte  is  er  een 
drie-dimensionaal  model  gebruikt,  waarbij  de  oneindige  hoogte  werd  vervangen  door  een  eindige 
(figuur  3). 

Samengevat: 

-  enkel  vat,  twee-dimensionaal  model,  oneindig  lange  cilinder, 

-  5x5  vaten,  twee-dimensionaal  model,  oneindig  lange  dlinders; 

-  idem  zonder  speciale  warmteoverdrachtselementen; 

-  2x5x5  vaten,  drie-dimensionaal  model,  eindige  lengte  van  dlinders. 
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3.3  Toegepaste  fysische  gegevens  van  her  kruit 

De  gegevens  uit  dit  onderzoc'’,  die  in  tabel  1  staan  vermeld,  zijn  afkomstig  uit  het  bestand  voor 
kruitgegevens  van  het  PML-TNO.  Het  kruit  met  aanduiding  "A"  is  eerst  toegepast.  Omdat  de 
resultaten  van  de  berekeningen  met  dit  kruit  niet  goed  te  verklaren  waren,  en  niet  meer  achter- 
haald  kon  worden  van  welke  soort  kruit  deze  gegevens  afkomstig  waren,  zijn  de  berekeningen 
herhaald  met  de  gegevens  van  kruit  "B".  Kruit  MC"  heeft  66n  van  de  hoogste  warmteprodukties  van 
de  kruiten  die  in  het  bestand  voorkomen.  Daarmee  is  het  meest  ongunsdge  geval  afgedekt. 

Tabel  1  Gegevens  over  de  kruiten 


Aanduiding  in  dit  rapport  van  kruitsoon. 

"A" 

"B” 

»c» 

Code  van  het  kruit 

(onbekend) 

KB  6954 

KB  5177 

Activeringsenergie  E 

[kj/mol] 

130,5 

122,7 

120* 

Warmteprod.  factor 

[W/kg] 

0,86xl021 

5,782xl016 

l,608x17 

Calorische  waarde 

[kj/kg] 

3037 

3037 

(onbekend) 

Soonelijke  warmte  Cp 

(J/kgK] 

1500 

1350* 

1350* 

Dichtheid  p 

[kg/m3] 

840 

840 

840* 

Warmtegeleidingsco6ffici6nt  1  [W/mK] 

0,08 

0,08 

0,08* 

Dit  gegeven  was  onbekend  voor  deze  soon  kruit;  daarom  is  een  kenmerkende 
waarde  voor  kruit  genomen. 


3.4  Gebruikte  hard-  en  software 

De  simulaties  zijn  uitgevoerd  met  het  eindige  elementen  pakket  Abaqus,  versie  4.9.  Dit  pakket 
draait  op  de  Convex  supercomputer  van  het  FEL-TNO,  Den  Haag. 

De  nauwkeurigheid  van  de  berekening  met  Abaqus  kan  worden  belnvloed  met  enkele  parameters. 
Voor  de  tempera tuur  is  0,1  Kelvin  de  opgegeven  tolerantie  (TEMTOL).  De  temperatuur  mocht 
op  elk  punt  maximaal  3  K  per  rekenstap  stijgen  (DELTMX). 
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4  RESULTATEN 

4. 1  Invloed  van  stapeling  op  de  kritische  omgevingstemperatuur 

De  kritische  omgevingstemperatuur  is  gedefinieerd  als  de  hoogste  omgevingstemperatuur,  waarbij 
een  hoeveelheid  kruit  voor  onbeperkte  tijd  kan  worden  opgeslagen  zonder  een  wegloopreactie  te 
vertonen.  Deze  grootheid  wordt  in  dit  onderzoek  gebruikt  als  maat  voor  de  stabiliteit  van  een 
bepaald  soort  kruit  die  op  een  bepaalde  manier  is  opgeslagen.  Daarbij  gaat  het  om  de  stabiliteit 
van  een  stapel  kruitvaten  in  vergelijking  met  66n  enkel  vat  kruit. 

Daartoe  is  eerst  66n  enkel  vat  gemodelleerd  (alle  afrnetingen  en  fysische  gegevens  staan  in  hoofd- 
stuk  3  beschreven).  Bij  een  dergelijke  eenvoudige  geometrie  kan  de  kritische  omgevings- 
temperar  ur  ook  analytisch  berekend  worden  (zie  hoofdstuk  2).  Ter  controle  van  het  model  is  de 
waarde  die  met  behulp  van  simulaties  is  bepaald  vergeleken  met  deze  analytisch  berekende 
waarde.  Het  resultaat  hiervan  staat  in  tabel  2.  Hieruit  blijkt  dat  de  met  simulaties  bepaalde 
kritische  temperatuur  ongeveer  een  graad  te  hoog  is.  Die  nauwkeurigheid  is  voldoende  om  bruik- 
baar  te  zijn  voor  het  vergelijken  van  verschillende  opslagconfiguraties  en  voor  het  beoordelen  van 
de  veiligheid  van  opslagsituaties. 

Tabel  2  Vergelijking  van  simulatie  (numeriek)  en  theorie  (analytisch) 

Kritische  omgevingstemperatuur  in  [*CJ  voor  een  enkel  vat  kruit* 


kruitsoort: 

"A" 

"B" 

"C" 

theorie 

30,8 

81,0 

64,3 

simulade 

31,8 

82,0 

65,0 

Voor  de  kenmerken  van  de  verschillende  soorten  kruit  wordt  verwezen  naar 
tabel  1  in  hoofdstuk  3 

De  kritische  omgevingstemperatuur  voor  kruittype  "A"  is  onrealistisch  laag.  Waarschijnlijk  zijn  de 
gebruikte  fysische  gegevens  niet  juist.  Daarom  is  de  berekening  herhaald  voor  twee  andere  kruit¬ 
type  n.  Kruittype  "B"  is  willekeurig  gekozen  uit  de  kruit-databank  van  het  PMl^TNO.  Om  het 
meest  ongunstige  geval  te  kunnen  onderzoeken,  is  het  kruit  met  een  van  de  hoogste  warmte- 
produkdes  uit  deze  databank  gekozen  als  kruittype  "C". 


Vervolgens  is  het  resultaat  van  de  simulade  van  een  enkel  vat  vergeleken  met  de  simulatie  van  een 
stapel  vaten.  Als  tussenstap  is  een  twee-dimensionaal  model  van  vijf  bij  vijf  vaten  gebruikt.  Dit 
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komt  overeen  met  een  oneindig  hoge  stapel.  Tenslotte  is  een  drie-dimensionaal  model  van  twee 
lagen  met  vijf  bij  vijf  vaten  doorgerekend.  Daarbij  is  aangenomen  dat  deze  lagen  zonder  tussen- 
ruimte  op  elkaar  staan.  Voor  de  lucht  tussen  de  vaten  is  een  warmtegeleidingscoSfficifnt 
van  1.0  W/mK  aangenomen.  Het  resultaat  van  al  deze  simulades  staat  in  tabel  3. 

Tabel  3  Kridsche  temperaturen  bepaald  door  simulaue 
Kriusche  omgevingstemperatuur  in  [*C] 


,.A« 

kruitsoort: 

"B" 

"C" 

aanta)  vaten 

1 

31,8 

82,0 

65,0 

5x5 

19,3 

63,9 

48,5 

2x5x5 

22,5 

68,6 

53,1 

Hoe  gevoeliger  het  opgeslagen  kruit  is  voor  thermische  invloeden,  des  te  lager  is  de  kridsche 
omgevingstemperatuur.  Zoals  reeds  gesteld,  komen  de  resultaten  van  kruitsoort  "A"  niet  overeen 
met  een  realisdsche  kruitsoort,  maar  ze  zijn  wel  geschikt  om  de  verschillende  modellen  met  eL  aar 
te  vergelijken.  De  kridsche  omgevingstemperatuur  van  een  stapel  met  een  zeer  gevoelig  kruit,  "C", 
komt  uit  rond  de  50  *C.  Bij  de  modellen  van  1  en  van  vijf  bij  vijf  vaten  vindt  de  afkoelL.g 
uitsluitend  plaats  via  de  zijwanden  (dit  komt  overeen  met  oneindig  hoge  stapels).  Door  vijf  bij  vi)'- 
vaten  tegen  elkaar  te  zetten  daalt  de  kridsche  omgevingstemperatuur  zo'n  zeventien  graden  ten 
opzichte  van  het  enkele  vat.  Het  drie-dimensionale  model  van  een  stapel  met  twee  lagen  is  wat 
realistischer,  omdat  hier  ook  afkoeling  via  de  boven-  en  onderwanden  plaats  heeft.  Dat  is  in 
werkelijkheid  ook  zo.  Daardoor  wordt  de  kridsche  omgevingstemperatuur  circa  vijf  graden  hoger 
dan  bij  de  oneindig  hoge  stapel  (5x5  model) .  In  werkelijkheid  is  de  afkoeling  zelfs  nog  iets  beter, 
omdat  er  onder  elke  laag  een  pallet  zit,  waarlangs  lucht  kan  stromen.  In  de  prakdjk  is  het  ook 
mogelijk  dat  de  vaten  een  klein  stukje  van  elkaar  af  staan.  Door  deze  verschijnselen  zal  dan  de 
kridsche  omgevings-temperatuur  nog  wat  hoger  uitkomen. 
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In  tabel  4  zijn  de  kritische  tempera  ruren  van  enkel  vat  en  een  stapel  met  2x5x5  vaten  nog- 

maals  onder  elkaar  gezet  om  het  verschil  te  laten  zien.  Hieruit  blijkt  dat  de  kritische  omgevings- 
temperatuur  bij  de  hier  gesimuleerde  omstandigheden  tot  circa  veertien  graden  lager  ligt  dan  de 
kritische  omgevingstemperatuur  voor  een  enkel  vat. 

Tabel  4  Invloed  van  stapeling  op  de  kritische  omgevingstemperatuur 


configurade: 

"A" 

kruitsoort: 

"B" 

"C" 

enkel  vat 

(Tx): 

31,8 

82,0 

65,0 

stapel 

(Ts): 

22.5 

68.6 

52A 

verschil 

(Ti-Ts): 

9,3 

13,4 

11.9 

Om  de  warmtehuishouding  in  een  stapel  kruitvaten  nader  te  beschouwen  zijn  de  tempera  tuur- 
profielen  van  een  simulatie  grafisch  weergegeven  in  figuur  4.  Het  gaat  hier  om  het  uiteindelijke  sta- 
biele  temperatuurprofiel  van  kruit  "C"  in  een  stapel  van  2x5x5  vaten,  bij  een  omgevings¬ 
temperatuur  van  53*C.  Deze  temperatuur  ligt  vlak  onder  de  kritische  omgevingstemperatuur 
(53, l'C).  Duidelijk  is  te  zien,  dat  het  middelste  punt  van  de  stapel  het  heetste  is  (dit  punt  ligt  aan 
de  rand  van  het  model,  omdat  slechts  een  representatief  segment  van  de  stapel  gesimuleerd  is).  De 
met  lucht  gevulde  kolommen  tussen  de  vaten  zorgen  voor  een  vergrote  warmteafvoer  naar  het 
bovenvlak  van  de  stapel. 

In  figuur  5  is  het  vetioop  van  de  temperatuur  tegen  de  djd  uitgezet,  voor  dezelfile  simulatie  als  bij 
figuur  4  en  tevens  voor  drie  simulaties  die  een  wegloopreactie  te  zien  geven.  Aangegeven  is  de  tem¬ 
peratuur  in  de  heetste  positie  in  de  stapel:  het  centrum  van  het  middelste  vat. 

In  alle  gevallen  loopt  de  temperatuur  de  eerste  den  dagen  circa  5  graden  op.  Bij  de  simulatie  met 
een  omgevingstemperatuur  van  55*C  loopt  de  temperatuur  daama  snel  op,  zodat  er  binnen  de 
volgende  tien  dagen  een  runaway-reactie  ontstaat.  Bij  de  simulaties  met  een  omgevings¬ 
temperatuur  van  54  en  53,5*C  neemt  de  tempera tuurstijgsnelheid  eerst  nog  wat  af,  maar  na  een 
aantal  maanden  treedt  er  toch  nog  een  runaway  op.  De  vierde  curve  hoort  bij  een  omgevings¬ 
temperatuur  van  53’C.  Dit  is  beneden  de  kridsche  omgevingstemperatuur,  waardoor  de  situatie 
stabiel  blijft. 
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TEmp 

VALUE 

1  +5.40E+01 

2  +5.50E+01 

3  +5.60E+01 

4  +5.70E+01 

5  +5.B0E+01 

6  +5.90E+01 

7  +6.00E+01 

0  +6.10E+01 

9  +6.20E+01 

10  +6.30E+01 

11  +6.40E+01 

12  +S.S0E+0 1 

13  +6.60E+01 


Figuur  4b  Contoutplot  van  een  stabiei  tempcraraurprofiel  (in  horizontaal  vlak  halvcrwege  de 
bovenste  laag  vaten) 


Figuur  4c  Contourplot  van  een  stabiei  tempera tuurprofiel  (in  verticaal  vlak  midden  door  het 
midden  van  de  stapel) 
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Figuur  5  Temperatuurverloop  bij  verschillende  omgevingstemperaturen  (model:  2x5x5 
vaten,  kruit:  KD  5177) 


4.2  Verfijning  van  de  berekening  van  de  kritische  temperatuur 

De  kritische  omgevingstemperaruur  van  een  opgeslagen  hoeveelheid  kruit  kan  als  volgt  bepaald 
worden.  Er  wordt  een  reeks  simulaties  uitgevoerd  met  verschillende  omgevingstemperaturen.  Als 
er  na  een  zekere  simulatieperiode  een  stabiel  tempera tuurprofiel  blijkt  te  ontstaan,  moet  er  een 
hogere  omgevingstemperaruur  geprobeerd  worden.  Als  er  daarentegen  een  runaway  optreedt, 
moet  de  omgevingstemperaruur  lager  gekozen  worden.  Deze  methode  is  erg  tijdrovend,  als  de  kri¬ 
tische  omgevingstemperatuur  nauwkeurig  berekend  moet  worden. 

Om  de  bepaling  te  versnellen,  kan  gebruik  gemaakt  worden  van  de  berekende  tijd  tot  runaway 
(inductietijd  tj).  Het  is  namelijk  gebleken,  dat  de  inverse  van  de  inductietijd  (1/tj)  min  of  meer 
lineair  varieert  met  de  omgevingstemperatuur  (zie  figuur  6).  Door  (1/tj)  uit  te  zetten  in  een  grafiek 
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en  te  extrapoleren  naar  ( 1/tj  =  0)  kan  een  schatting  van  de  kritische  omgevingstemperatuur  worden 
gevonden.  Deze  kan  als  leidraad  dienen  voor  de  keuze  van  de  omgevingstemperatuur  voor  de 
volgende  simulade.  Daardoor  kan  het  aantal  benodigde  iterades  beperkt  worden. 


7.00£a)7 

6.00647 

$00647 

_  4  00647 

6 

”  3  00647 

200647 

100647 

0.006*00 


4$  46  47  46  40  SO  51  52  S3  $4 

T»PO 


Figuur  6  Relade  tussen  de  inducdeujd  en  de  omgevingstemperatuur 


4.3  Gevoeligheidsanalyse 

In  deze  paragraaf  wordt  nagegaan  wat  de  invloed  van  verschillende  parameters  is  op  de  simulatie- 
resultaten.  Daarmee  wordt  een  idee  gekregen  van  de  nauwkeurigheid  van  de  resultaten.  Ook  wordt 
daarmee  duidelijk  hoe  gevoelig  het  resultaat  is  voor  een  variade  in  de  waarde  van  verschillende 
parameters. 
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4.3.1  Effect  van  het  gekozen  rekenrooster 

De  invloed  van  het  gekozen  rekenrooster  is  onderzocht  aan  de  hand  van  het  model  van  een  enkel 
vat.  Hiervoor  zijn  drie  verschillende  grids  gebruikt,  met  respectievelijk  adit,  twaalf  en  tweeindertig 
elementen  (zie  flguur  1).  De  gevonden  kritische  omgevingstemperatuur  staat  in  tabel  5.  Bij  deze 
berekening  zijn  fysische  parameters  gebruikt  van  een  denkbeeldige  kruitsoort,  waardoor  de 
gevonden  kritische  tempera turen  als  zodanig  niet  representatief  zijn.  Uit  het  resultaat  blijkt,  dat 
gridverfijning  niet  nodig  is.  Voor  het  model  van  een  enkele  drum  is  verder  het  grid  van 
tweefindertig  elementen  gebruikt.  Het  grid  voor  de  modellen  van  een  stapel  vaten  is  wel  gebaseerd 
op  het  grofste  grid,  teneinde  het  aantal  rekenpunten  niet  te  groot  te  laten  worden. 

Tabel  5  Kritische  omgevingstemperatuur  bij  verschillende  rekenroosters 


Kritische  omgevingstemperatuur  [*C] 

grof  grid 

(8  elementen) 

22.4 

fijngrid 

(12  elementen) 

22.4 

extra  fijn  grid 

(32  elementen.) 

22.6 

4.3.2  Invloed  van  de  warmtegeleidingsco£ffirient  van  kruit 

Met  behulp  van  de  thermische  explosie  theorie  (zie  hoofdstuk  2)  kan  de  kritische  omgevings¬ 
temperatuur  voor  een  oneindig  lange  cilinder  op  eenvoudige  wijze  worden  berekend.  Daarmee  kan 
ook  de  invloed  van  de  warmtegeleidingscoCfficient  van  kruit  onderzocht  worden.  Het  resultaat  van 
deze  berekeningen  staat  in  tabel  6.  De  fysische  gegevens  van  kruittype  "A"  zijn  toegepast. 

Tabel  6  Invloed  van  de  warmtegeleidingscoifSdSnt  van  kruit 


warmtegeleidingsco€ffici£nt  kruit 

kritische  omgevingstemperatuur 

W/mK] 

rcj 

0,08 

30,8 

0,10 

32,2 

0,12 

33,3 

0,15 

34,7 
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Uit  deze  resultaten  blijkt  dat  de  warmtegeleidingscoeffidCnt  van  kruit  een  belangrijke  invloed  op 
de  kritische  omgevingstemperatuur  heeft.  Een  variatie  van  de  warmtegeleid  ingscofifficient  van 
0,08  tot  0,15  W/mK  gee  ft  circa  vier  graden  verschil  in  de  kritische  omgevingstemperatuur. 

4.3.3  Invloed  van  de  warmteoverdrachtsco£fEci£nt 

Er  vindt  warmteoverdracht  plaats  tussen  het  kruit  en  het  verpakkingsmateriaal  en  tevens  tussen 
verpakking  en  de  lucht  in  de  omgeving.  De  warmteflux  door  het  buitenoppervlak  is  evenredig  met 
het  temperatuurverschil  tussen  de  buitenrand  van  het  kruit  en  de  lucht  er  omheen.  De  even- 
redigheidsconstante  is  de  warmteoverdrachtscoCfficiCnt  en  voor  vrije  convectie  is  de  waarde 
daarvan  tussen  de  3  en  de  20  W/Km2  [3).  De  warmteoverdrachtscoCfficiint  heeft  geen  grote 
invloed  op  de  kritische  omgevingstemperatuur,  want  de  grootte  van  de  warmtestroom  wordt 
vooral  bepaald  door  de  warmtegeleidbaarheid  van  het  kruit  zelf.  In  dit  onderzoek  is  een  warmte- 
overdrachtscofifficiCnt  van  6  W/km2  genomen.  Om  de  invloed  hiervan  nader  te  onderzoeken,  zijn 
ook  enkele  simulates  uitgevoerd  met  een  twee  maal  zo  grote  warmteoverdrachtsco£ffici£nt.  De 
resultaten  hiervan  staan  in  tabel  7.  De  fysische  gegevens  van  kruittype  "A"  en  het  model  van  een 
enkel  kruitvat  zijn  gebruikt  (zie  hoofdstuk  3).  De  conclusie  is  dat  de  invloed  van  de  warmte- 
overdrachtscofifficient  zeer  gering  is. 

Tabel  7  Invloed  van  de  warmteoverdrachtscoCfficiCnt 


warmteoverdrachtscoCffiriCnt 

kritische  omgevingstemperatuur 

[W/m2K] 

[•ci 

6.0 

31.8 

12.0 

32.1 

4.3.4  Invloed  van  interface-elementen  in  het  model 

In  het  oorspronkelijk  model  van  5x5  kruitvaten  zijn  interface-elementen  opgenomen,  waarmee  de 
warmteoverdracht  tussen  het  kruit  en  de  lucht  tussen  de  kruitvaten  wordt  uitgerekend.  Hierboven 
is  reeds  vermeld,  dat  de  warmteweerstand  voomamelijk  bepaald  wordt  door  de  warmtegeleid¬ 
baarheid  van  het  kruit  zelf,  zodat  de  waarde  van  de  warmteoverdrachtscoCfficiCnt  weinig  invloed 
op  de  resulterende  warmtestromen  heeft.  Daarom  zijn  deze  interface-elementen  uit  het  model 
verwijderd,  ten  einde  de  benodigde  rekendjd  te  beperken.  De  invloed  van  deze  vereenvoudiging  op 
de  kritische  omgevingstemperatuur  is  te  zien  in  tabel  8.  De  gebruikte  kruitgegevens  zijn  van 
kruittype  "A"  en  het  model  is  van  5x5  kruitvaten  (met  en  zonder  interface-elementen,  zie 
hoofdstuk  3). 
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Tabel  8  Invloed  van  de  interface-el ementen  in  het  model 


model 

kritische  omgevingstemperatuur  [*C] 

met  interface-elementen 

19,2 

zonder  interface-elementen 

19,3 

De  interface-elementen  kunnen  dus  zonder  noemenswaardig  verlies  aan  nauwkeurigheid  uit  het 
model  weggelaten  worden.  De  resultaten  van  paragraaf  4.1  zijn  verkregen  met  het  model  zonder 
interface-elementen.  De  simulades  die  in  paragraaf  4.3.5  worden  beschreven  waren  reeds  uit- 
gevoerd  met  het  model  mfct  interface-elementen.  In  het  drie-dimensionale  model  zijn  geen 
interface-elementen  toegepast. 

4.3.5  Invloed  effectieve  warmtegeleidingscoeffiriSnt  van  lucht 

In  hoofdstuk  2  is  uiteengezet,  dat  de  warmtestromen  door  de  luchtruimtes  tussen  de  kruitvaten 
met  de  hier  gebruikte  methode  niet  exact  te  beschrijven  zijn.  Er  is  daarom  een  benaderende  me- 
thode  gebruikt,  waarbij  aan  de  lucht  een  veihoogde  warmtegeleidingscofifficiCnt  is  toegekend. 
Daarbij  is  het  warmtetransport  ten  gevolge  van  de  beweging  van  de  lucht  mede  vendisconteerd  in 
de  warmtegeleidingscofiflBd5nt.  De  warmtegeleidingscoefficiCnt  van  een  sugnante  (onbeweeglijke) 
luchtlaag  is  0,024  W/mK  (bij  273  K)(3].  De  wanmegeleidingscoCffidCnt  van  kruit  ligt  ongeveer 
tussen  0,08  en  0,15  W/mK.  In  het  model  waarmee  de  overige  berekeningen  zijn  gedaan,  is  de 
effectieve  warmtegeleidingsco£fSci£nt  van  lucht  op  1,0  W/mK  gesteld.  Om  de  .nvloed  van  deze 
parameter  te  onderzoeken,  zijn  een  aantal  simulades  met  een  hogere  of  lagere  effectieve 
wanntegeltidingsco£fBci£nt  van  lucht  uitgevoerd.  De  resultaten  hiervan,  die  zijn  verkregen  voor 
kruittype  "A"  en  het  model  met  5x5  kruitvaten  (zie  hoofdstuk  3),  staan  in  tabel  9.  Ter  vergelijking: 
de  kritische  omgevingstemperatuur  voor  een  enkel  vat  komt  op  31,8*C. 
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Tabei  9  Invloed  van  de  effectieve  warmtegeleidingsco€ffici£nt  van  lucht 


effectieve  warmte-geleidingsco££fici£nt 
[W/mK] 

kritische  omgevingstemperatuur 

rcj 

0,08 

15,5 

1,0 

19,2 

10,0 

22,3 

Uit  deze  resultaten  kan  geconcludeerd  worden  dat  de  effectieve  warmtegeleidingsco€fBci£nt  een 
belangrijke  invloed  heeft  op  de  kritische  omgevingstemperatuur.  Als  een  tien  maal  zo  lage  waarde 
voor  de  wanntegeleidingscoefficient  van  lucht  wordt  genomen,  komt  de  kritische  omgevings- 
temperatuur  circa  vier  graden  lager  uit. 

Een  overzicht  van  de  uitgevoerde  simulaties  staat  in  bijlage  1. 

5  CONCLUSffiS  EN  AANBEVEUNGEN 


5.1  Conclusies 

Bij  de  opslag  van  kruit  kan  er  ontbranding  ten  gevolge  van  zelf-opwanning  plaatsvinden,  indien  de 
omgevingstemperatuur  boven  een  kritische  waarde  ligt.  Deze  kritische  omgevingstemperatuur 
hangt  niet  alleen  van  de  eigenschappen  van  het  kruit  af,  maar  ook  van  de  omvang  van  de 
verpakkingseenheden  en  de  wijze  waarop  deze  gestapeld  zijn.  In  het  in  dit  rapport  beschreven 
onderzoek  is  de  invloed  van  stapeling  op  de  kritische  omgevingstemperatuur  onderzocht  door 
middel  van  numerieke  Simula  tie.  Door  middel  van  deze  simulaties  is  een  beeld  gevormd  van  de 
warmtestromen  in  een  opgeslagen  hoeveelheid  kruit.  Het  is  t evens  duidelijker  geworden  hoe  het 
tempera tuumiveau  waarbi)  kruit  veilig  kan  worden  opgeslagen  belnvloed  wordt  door  de  wijze  van 
opslag. 
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De  nu  volgende  conclusies  zijn  van  toepassing  op  enkele  kruitsoonen  die  in  een  veel  gebruikt  type 

cilindervormige  kruitvaten  opgeslagen  zijn,  op  een  stapel  van  twee  lagen  van  vijf  bij  vijf  vaten,  zoals 

beschreven  in  hoofdstuk  3. 

-  Bij  de  onderzochte  stapel  kruitvaten  bleek  de  kritische  omgevingstemperatuur  tien  tot  vijftien 
graden  lager  dan  de  kritische  omgevingstemperatuur  van  een  enkel  vat. 

-  Voor  een  kruitsoort  met  een  relatief  hoge  warmteproduktie  (KD  5177)  ligt  de  kritische 
omgevingstemperatuur  voor  een  stapel  vaten  met  kruit  rond  de  50*C.  Bij  een  andere 
onderzochte  kruitsoon  (KD  6954)  ligt  deze  waarde  rond  de  65  *C. 

-  De  waarde  van  de  effectieve  warmtegeleidingscofcffitient  van  de  ludit  tussen  de  vaten  heeft  een 
grote  invloed  op  de  kritische  omgevingstemperatuur.  Bij  de  meeste  simulaties  is  een  waarde 
van  1,0  W/mK  aangenomen.  Bij  een  waarde  van  0,08  W/raK  werd  de  kritische  omgevings¬ 
temperatuur  circa  vier  graden  verlaagd.  Bij  een  waarde  van  10  W/mK  werd  de  kritische 
omgevingstemperatuur  daarentegen  drie  graden  hoger. 

-  De  inverse  van  de  inductietijd  of  ujd-tot-runaway  bij  temperaturen  die  iets  boven  de  kritische 
omgevingstemperatuur  liggen  blijkt  min  of  meer  lineair  te  variCren  met  de  omgevings¬ 
temperatuur.  Dit  verschijnsel  kan  toegepast  worden  om  het  aantal  benodigde  iteraties  bij  het 
numeriek  bepalen  van  de  kritische  omgevingstemperatuur  te  beperken. 

-  De  warmtegeleidingscoCfficiCnt  van  kruit  heeft  een  belangrijke  invloed  op  de  kritische 
omgevingstemperatuur.  Een  variatie  van  de  warnuegeleidingscoifficiCm  van  0,08  tot 
0,15  W/mK  geeft  circa  vier  graden  verschil  in  de  kritische  omgevingstemperatuur. 

-  De  grootte  van  de  warmteoverdrachtsco£ffici£nt  voor  de  warmteoverdracht  tussen  het  vat  en  de 
omgeving  heeft  weinig  invloed  op  de  kritische  omgevingstemperatuur.  Er  zijn  twee  simulaties 
uitgevoerd  met  een  warmteoverdrachtscoCfficiCnt  van  respectievelijk  6  en  12  W/m2K.  Het 
verschil  in  kritische  omgevingstemperatuur  dat  daaruit  voortvloeide  was  slechts  0,3’C. 
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5.2  Aanbevelingen 

-  In  het  hier  beschreven  onderzoek  is  de  temperatuur  in  de  omgeving  van  de  stapel  als  een 
constante  beschouwd.  De  totale  voorraad  kruit  kan  echter  een  aanzienlijke  hoeveelheid  warmte 
afgeven,  waardoor  de  omgevingstemperatuur  zou  kunnen  stijgen.  Om  dit  na  te  gaan  is  het 
nutdg  om  de  maximale  warmteprodukde  van  de  stapel  uit  te  rekenen  bij  een  omgevings- 
temperatuur  in  de  buurt  van  de  kridsche  waarde. 

-  Het  is  nutdg  om  na  te  gaan  of  het  gedrag  van  een  stapel  kruitvaten  ook  met  een  veel 
eenvoudiger  model  ingeschat  kan  worden.  Daarvoor  zou  de  stapel  als  een  enkel  homogeen  blok 
gemodelleerd  moeten  worden-  waarbij  de  fysische  parameters  op  een  bepaalde  manier  uit  de 
parameters  van  kruit  en  lucht  afgeleid  worden.  Hiermee  zou  ook  weer  de  kridsche  omgevings- 
teuperatuur  en  de  tijd-tot-runaway  voor  verschillende  omstandigheden  uitgerekend  kunnen 
worden. 

-  Tot  nu  toe  is  aangenomen  dat  de  begintemperatuur  van  het  kruit  gelijk  is  aan  de  omgevings¬ 
temperatuur.  Dat  hoeft  in  de  prakdjk  niet  aldjd  het  geval  te  zijn,  bijvoorbeeld  als  het  kruit  nog 
warm  is  van  het  produkdeproces.  Zo  blijkt  uit  een  zeer  recente  publicade  [4],  dat  het  plaatsen 
van  een  warm  vat  eneigedsch  materiaal  in  een,  ter  compensade,  koele  cel  maar  een  zeer  gering 
effect  te  hebben  op  de  kridsche  temperatuur  van  het  vat:  ongeveer  1  l’C  daling  van  de  kridsche 
temperatuur  bij  een  temperatuurdaling  van  140*C  van  de  omgeving!  Deze  berekeningen 
werden  uitgevoerd  voor  een  enkel  kubisch  vat.  Aanbevolen  wordt  om  de  invloed  stapeling  op 
interne  kridsche  vattemperatuur  in  relatie  met  de  omgevingstemperatuur  door  middel  van 
simulaues  nader  te  onderzoek  en. 

-  Om  de  warmtestromen  in  de  luchdagen  in  de  stapd  nog  beter  te  onderzoeken  is  experimented 
onderzoek  wenselijk.  Dit  onderzoek  zou  ondermeer  moeten  resulteren  in  een  betrouwbare 
schatting  van  de  effectieve  warmtegelddingscoCfficiCnt  van  lucht  ten  behoeve  van  simulaties. 
Als  het  experiment  vooral  op  het  onderzoeken  van  de  warmtestromen  tussen  de  vaten  is 
gericht,  kan  het  uitgevoerd  worden  met  een  inerte  stof  voorzien  van  dektrische  verw arming. 
Daarmee  worden  problemen  met  de  veiligheid  voorkomen. 
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BIJIAGE  1  OVERZICHT  VAN  SIMULATIERUNS 

In  de  volgende  tabel  zijn  van  de  uitgevoerde  simulatieruns  de  volgende  gegevens  opgenomen: 

-  naam  van  de  simulatie-run  (volgnummer) 

-  model : 

la  :  model  van  1  vat  (bestaande  uit  8  elementen) 

lb  :  model  van  1  vat  (bestaande  uit  12  elementen) 

lc  :  model  van  1  vat  (bestaande  uit  32  elementen) 

55  :  model  van  5  bij  5  vaten 

55z:  model  van  5  bij  5  vaten  zonder  interface-elementen 

255:  model  van  2  lagen  van  5  bij  5  vaten 

-  kruitsoort: 

"A",  MB”  of  "C"  :  zie  hoofdstuk  3, 

"X":  acdveringsenergie  115  kj/mol; 

warmteproduktiefactor  0,775x10^  W/kg; 
dichtheid:  840  kg/m3; 
warmtegeleidingsco£f5ci£nt:  0,08  W/mK; 
soortelijke  warmte:  1300  J/kgK. 

-  omgevingstemperatuur  (T^) 

-  optreden  runaway  (r.a.): 

+:  wel  een  runaway; 

- :  geen  runaway; 

+/-:  twijfelgeval. 

-  maximaal  bereikte  tempera tuur  in  warmste  punt  (Tmax) 

-  tijd-tot-runaway  (inductietijd)  in  uren  (tj) 

-  opmerkingen  : 

Tj:  ini  tide  tempera  tuur  van  het  kruit 

h:  warmteoverdrachtsco&ffitiCnt 

1):  effectieve  warmtegelekUngscottEciCnt  van  lucht 
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Tabel  1.1  Overzicht  van  simuladeruns  (vervolg) 


run 

model 

kruit 

Ta  [*C] 

r.a 

Tmax  t'C] 

t|  uur 

opmerkingen 

kalO 

55 

"A" 

15 

- 

20,6 

1,  =  0,08  W/mK 

ka9 

55 

"A" 

16 

+ 

6960 

1]  =  0,08  W/mK 

ka7 

55 

"A" 

18 

+ 

1368 

ll  =  0,08  W/mK 

kcl 

55z 

"A" 

19 

- 

26,3 

kc2 

55z 

"A" 

20 

+ 

2717 

kc3 

55z 

"A" 

21 

+ 

1146 

kd3 

255 

"A" 

21 

- 

25,0 

kd7 

255 

"A" 

22,5 

- 

30,0 

kd8 

255 

"A" 

23 

+ 

1556 

kd6 

255 

"A" 

24 

+ 

530 

kd5 

255 

"A" 

25 

+ 

333 

kc6 

lc 

"B" 

82 

- 

? 

kc4 

lc 

-B" 

82,5 

- 

±583 

kc2 

lc 

„B» 

84 

- 

140 

ke3 

ic 

"B" 

86 

- 

72,8 

kf2 

55z 

"B" 

62 

- 

67,8 

kf5 

55z 

"B" 

63,5 

- 

73,8 

kf4 

55z 

„B« 

64 

+/- 

79,0 

kf3 

55z 

"B" 

65 

+ 

2314 

kfl 

55z 

"B" 

66 

+ 

1205 

kgl 

255 

"B“ 

68 

- 

75,7 

kg3 

255 

"B" 

69 

+/- 

81,6 

k«2 

255 

"B" 

70 

+ 

875 

kg4 

255 

«B« 

71 

+ 

505 

kh5 

lc 

"C" 

64,5 

- 

73,5 

kh4 

lc 

"C" 

65 

+/- 

77,3 

khl 

lc 

"C" 

66 

+ 

287 

kh2 

lc 

"C" 

67 

+ 

140 

kh3 

lc 

"C" 

68 

+ 

93,4 
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Tabel  1.1  Overzicht  van  simuladeruns  (vervolg) 


run 

model 

kruit 

Ta  [*C1 

Tmaxl'QI 

tj  uur 

opmcrkingen 

ki7 

55z 

"C" 

43,5 

- 

46,3 

ki5 

55z 

"C" 

48 

- 

61,4 

ki4 

55z 

"C" 

48,5 

+/- 

53083 

ki3 

55z 

"C" 

49 

+ 

4564 

kil 

55z 

*C" 

50 

+ 

1476 

ki2 

55z 

»C» 

51 

+ 

928 

ki6 

55z 

"C" 

53,5 

+ 

431 

kkl 

255 

"C" 

53 

- 

65,0 

kk4 

255 

"C" 

53,5 

+ 

1886 

kk2 

255 

"C" 

54 

+ 

905 

kk3 

255 

"C" 

55 

+ 

491 
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